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ABSTRACT: Hadron multiplicities — especially for strange particles — are calcula- 
ted in the framework of the algebraic coalescence rehadronization model (ALCOR), 
which counts for redistribution of quarks into hadrons for relativistic heavy-ion 
collisions. The influence of Bjorken flow on the final hadronic composition are 
incorporated in the model. A comparison is made with the CERN SPS NA35 S + S 
and p + p experiments. The analysis of these experiments with ALCOR shows a 
strangeness enhancement for S-|-S collisions and a possible formation of a sort of 
semi-deconfined state of the matter. Predictions for Pb-|-Pb collisions (NA49) are 
also presented. 



Introduction 

The strongly dynamical process of rehadronization from a thermal quark-gluon 
plasma or from any other type of prehadronic quark matter into the finally 
observed hadrons in a high-energy heavy-ion collision cannot be described in 
the framework of perturbative parton dynamics alone [Q] , [0] nor by using rate 
equations as in quarko- and hadrochemistry P],^. Not only does the strength 
of the interaction leading to valence quark confinement into hadrons in the final 
state of this process rule out any expansion in terms of the QCD coupling 
constant g'^/An, but also the fact that there are only color singlet hadrons 
renders a description using rate equations derived in the infinite volume limit 
to be unreliable. In this case there would always be a small fraction of free 
quarks present after a finite time has elapsed unless the hadronization cross 
section diverges after a finite time. In that case the use of rate equations is 
also unjustified. 

The idea of redistributing all the quarks and antiquarks present in the 
quark matter at the latest stage of its evolution between color singlet final 

^ The talk was presented by J. Zimanyi on S'95 Conference, Tucson, Jan. 3-7 1995. 
E-mail:jzimanyi@sunserv. kfki.hu 
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states (hadrons) cures these two problems at the level of a phenomenological 
model. It describes the rehadronization as a sudden process. This assumption 
unfortunately cannot be checked in the framework of a quark number redis- 
tribution model. For this purpose a comparison with microscopic dynamical 
calculations will be necessary. If one then compares the predictions of such 
a simple redistribution of different quark flavors with experimental findings, 
especially with strange or charm to light hadron ratios, one may then have a 
basis for judging whether nontrivial collective phenomena at the quark level 
must be considered in these experiments. 

In applying rehadronization models one has to supply as input the num- 
ber of light and heavy valence quarks and antiquarks, all of which will be 
hadronized. This number can be obtained from theoretical models (e.g. string 
model 0, color rope model [0, thermodynamical models, etc.) or from a part 
of the observed data. In the latter case comparison of the predictions of the 
model with further experimental data may be used to prove or disprove the 
basic assumptions of that model. Here, applying the ALgebraic COalescence 
Rehadronization model, ALCOR we shall follow the latter possibility. 

The ALCOR model 

The starting point of the ALCOR model is a state in which only those quarks 
and antiquarks which form the final hadrons are present in a local thermal equi- 
librium. We also assume that they are all present, so their number already 
account for some of the gluons fragmented earlier. This assumption is easily 
acceptable for the heavy quarks for which any meeting with their antiparticle is 
improbable due to their small specific density. In the case of light constituent 
quarks (counted after gluon fragmentation) the approximate entropy conser- 
vation during the rehadronization process leads to the conservation of their 
numbers. This picture was first formulated in the framework of the algebraic 
recombination model [§. In this model the number of a given type of hadron 
produced is proportional to the product of the numbers of their constituting 
quarks. This picture was developed further by the introduction of gluon frag- 
mentation 0. We shall denote the number of produced u,d,s quark pairs, 
just before the hadronization, by N^^pair, N^^pair and Ng^pair, respectively. The 
strangeness production factor, is defined as gs = Ng^pair/iNu^pair + N^^pair)- 
(The definition of factor fs in Ref. [0 was different from that of Qs-) 

In the ALCOR model we generalize this further by introducing the coa- 
lescence factors CM^hj), CB{i,j, k) for each hadron type separately. In order 
to ensure that all quarks are included in hadrons in the final state we also 
introduce normalization factors h{i) associated with each quark flavor when 
counting the quarks and antiquarks in hadrons. Furthermore, a spin degen- 
eracy factor, Z}(^) = 2Sh + 1, is also included. Here we shall consider the 
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simplest hadron multiplets conserving isospin symmetry. These are the spin 
and spin 1 meson octets together with the a and uj isospin singlet meson 
states. The spin 1/2 octet and spin 3/2 decouplet baryons and anti-baryons 
are also included. Thus the numbers of a given sort of baryon or antibaryon, 
consisting of quarks with flavors j, fc, are given as 

= D('^CB{t,J,k)b{t)b{j)b{k)NQ{^NQ{j)NQ{k), (1) 
Nf{z,j,k) = D^^'^C^{t,j,k)h{i)h{j)h{k)NQ{t)NQ{j)NQ{k), (2) 

and the number of a given sort of meson is given as 

Nthh3)=D^^^CM{^.3)h{mJ)NQ{i)NQ{3). (3) 

Here Nq{i) and N-q{i) are the number of quarks and antiquarks of type i , 
respectively. 

These constraints lead finally to the following system of equations (one for 
each quark and each antiquark flavor) 

N; 

h j=l 

N; Nf 

+ E E E(l + + kk)D^'^CB{z,J, k)bmj)b{k)NQ{z)NQ{j)NQ{k), 

h j=l k=j 

(4) 

Nf 

h j=l 

Nf Nf 

+ E E E(i + + s.,k)D^'^c^{z,j, kmmmNQmQmQik). 

h j=l k=j 

(5) 

The number of independent equations in eqs.(|D,(|D is equal to the number 
of independent b{i) and b{i) factors. Note, that these 2Nf quantities are 
not adjustable parameters of the model, but they are normalization constants 
determined by eqs.(4-5). 

The Coalescence Factors 

Let us consider now two particles (quarks or antiquarks) of type (flavor) 'a' 
and 'b' which form a particle (hadron) 'h'. The rate equation describing this 
formation of particle 'h' as a one-way process is based on the number density 
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of particles 'a' and 'b' {na,nb) in the reaction volume V and on the specific 
cross section a of the a + b ^ h process: 

^^i'^'n^h) = {crv)nanb . (6) 

Here {av) is the thermally averaged rate including the relative velocity v of 
the reacting partners. Starting such a process with no 'h' particles at the 
beginning the rate of change is initially constant — (o" f )na(ro)n5(ro) — thus 
the number of 'h' particles grows linearly. After a characteristic time r we 
have 

Nh = Vn^, = ^^^NaNb = C{a, b)N,Nb (7) 

particles. In a sudden approximation we see the effect of different branching 
ratios due to competition between several possible reaction channels producing 
different types of particles (protohadrons). In this way a coalescence factor 

C(a,6) = Ml (8) 

can be defined for each hadron. These factors infiuence the relative distribution 
of strange quarks between strange mesons and singly or doubly, etc. strange 
baryons. 

Considering the ratios of different hadron abundancies in the final state the 
proper time (r) and volume {V) dependence cancels, therefore the validity of 
the sudden approximation can not be checked by a direct comparison with ex- 
perimental data. A long hadronization time, r, on the other hand, would lead 
to chemical equilibrium between the hadrons. In this case the relative abun- 
dances would depend only on final state parameters and not on the branching 
ratios based on specific hadronization channels. It must be a very improba- 
ble coincidence if these two different assumptions about the rehadronization 
dynamics would lead to similar strange/non-strange hadron or meson/baryon 
ratios. Thus an indirect experimental check on the sudden approximation is 
provided. 

We assume furthermore that the creation of a baryon from three quarks or 
an antibaryon from three antiquarks is a two-step process leading through an 
intermediate diquark formation in the corresponding color triplet or antitriplet 
state which forms the final baryon later on together with a third quark or an- 
tiquark. These diquarks must, however, be very unstable, short-lived clusters 
(spatial correlations), so they may decay before forming the baryon. We take 
into account this mechanism of baryon production by introducing the phe- 
nomenological baryon production parameter, gs- 
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Thus the coalescence factor of a three-quark hadron is calculated as 



CB{a,b,c) = gB^{CM{a,b)CM{[a + b],c) + 

+CM{a, c)CM{[a + c], h) + Cuih, c)Cm([& + c], a)} , (9) 

where (/b is assumed to be flavor blind and we average on all combinatoric 
possibility to form a baryon N^\a, b, c). Here the arguments a, b and c stand 
for any of the flavors u, d or s considered, respectively. 

Note that the general scheme displayed in eqs.(|I] - 0), namely that the 
number of formed hadrons is proportional to the product of the numbers of 
initials, Nh = C{a,b)NaNb, has also been conjectured in the framework of a 



relativistic coalescence model |]10|- In this latter case, however, the coalescence 



factor, C(a, 6), has been obtained in another way: only an initial and final state 
wave function overlap in a Gaussian wave packet approximation was taken into 
account instead of using a transition matrix element inherent in the definition 
of the cross section involved in the ALCOR model. 

There are two ingredients occurring in eq.® which have not yet been 
discussed: i) the reaction volume, ii) the hadronization rate. The first is 
mainly restricted by the presence of a flow preventing some candidate valence 
quarks 'a' and 'b' from meeting to form a hadron 'h' if they are sitting in fluid 
cells which move away from each other so fast that the local thermal motion 
cannot overcome this relative velocity barrier. 



In ref. [|rT| the space-time and momentum-space evolution of the reaction 
at SPS energies was investigated using the SPACER model. It was found 
that the momentum distribution of the particles can be approximated by a 
local thermal distribution superimposed on a scaling flow in a finite space- 
time rapidity interval. We perform the calculation of the average reaction rate 
for this case. 



Hadronization Rate in the Presence of Flow 

In the presence of a flow we consider the relativistic Jiittner distribution 

=e-^P-"(^) (10) 

with u^{x) being the local four-velocity of the flow. Let us restrict ourselves 
in the following to a scaling longitudinal flow |jl2|. In this case the thermal 
averaging of the rate can be interpreted only locally: any rate equation of type 
eq.(|^) has to be reinterpreted in terms of reacting components from different 
coordinate rapidity ranges r]a and rji, respectively 

Id f 

— —{VUh) = / dT]adT]b{<Jv)abnanb. (11) 
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Clearly the total volume of the fireball increases due to a longitudinal scaling 
expansion like V = rnR"^ Ar] with total coordinate rapidity extension At] 
and transverse radius R. Assuming a space-time rapidity plateau for a finite 
interval in the number distribution of the reacting and produced components 
we arrive at a modified rate equation 



= -^N,N, . (12) 



A = / d(r]a - rjb) {av) Q{{pa - Pb) ■ {xa - Xb)) (13) 



Here 



is the total rate of reactions between all possible coordinate rapidity cell pairs. 
The constraint {pa — pt) ■ {xa — Xb) > reduces to the requirement that the 
relative velocity vector is oriented opposite to the relative position vector in the 
center of mass system of the colliding components 'a' and 'b'. This constraint 
is necessary because only those particles collide which have a relative velocity 
pointing towards each other. 

The inclusion of this constraint after some algebraic manipulation leads to 
the general expression of the averaged rate 

%mlmlK2{(3ma)K2{(3mb) ^ ^ 

where \ab{,s) = [s — {nia + mb)'^] [s — {rria — mb)"^] and G{Py/s) is the thermal 
weight factor. In the presence of longitudinal scaling flow the latter can be 
written as an integral over relative coordinate rapidity f] = {rja — rjb) 



K,{l3^sdi\7^\)-KJp^sdi\7^\ + /?^sh|r7|) 

G{(3^s) = d7]^ ^ ^. 

•^-oo pXj^ (s) chlr^l shlr^l 

(15) 

Investigating the integrand of eq. ([T5|) one observes a finite width, St], in the 
relative coordinate rapidity rj. (This quantity, 61], is equal to the relative flow 
rapidity because of the Bjorken scaling assumption.) The finite width depends 
on temperature T, on the particle rest masses ma,mb, and on the considered 
energy scale ^/s in a complicated manner. 
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Elementary Cross Sections 



Finally, the most uncertain ingredient of the ALCOR hadronization model, 
the hadronization cross section will be discussed. Here, until the dynamical 
confinement mechanism in QCD has been explored, an analogy with the p + 
A ^ d+{A — l) nuclear rearrangement (pick-up) reaction leading to deuteron 



formation may give some hints |T3[. Considering a final state in a Coulomb- 
like potential — a simplified picture of mesons at intermediate temperatures 
— a fusing cross section of 

2 

a = 16m?Vip3 — ^ (16) 



can be derived from the above mentioned analogy. Here rrih is the rest mass of 
the meson, while a = 1/ {niabd) is the Bohr radius of the bound qq state in a 
V{r) = —a/r Coulomb potential with niab being the reduced mass of particles 
'a' and 'b'. In our calculation we took a = 0.46 and T = 200 MeV temperature 
in estimating hadronization cross sections at CERN SPS energy. The factor 
p = 0.3 fm occurring in eq.([l6D accounts for the medium infiuencing the hadron 
formation and was taken to be equal to the Debye screening length in quark- 
gluon plasma at the above temperature. Finally k occurring in eq.(|l^) is the 
magnitude of the relative momentum vector of particles 'a' and 'b' measured 

1 /2 

in their center of mass system k = A^^ {s)/2^/s. 



Results 



Well armed with the above theoretical considerations the redistribution of dif- 
ferent fiavor quarks and antiquarks into all possible hadrons can be calculated 
in the ALCOR model using practically only three parameters: i) the total 
number of quark- antiquark pairs, Ntot,pair = Nu^pair + Nd,pair + Ns^pair, which 
can be determined from the measured total charged multiplicity, ii) the param- 
eter qb controlling the baryon formation and iii) the strangeness production 
factor Qg. The dependence of the final results on the other parameters (e.g. 
T, a) within their physically acceptable interval is small. All further results, 
such as the number of hyperons, kaons, etc. are predictions of the ALCOR 
model. 
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Table 1: Hadron multiplicities for S + S collision observed experimentally 
and obtained in the ALCOR and in the HIJING, RQMD, QGSM models 
(early non-collective string model versions, taken from refs. |[T^ |T5[) at 200 
GeV/nucleon bombarding energy. 
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Table 1 shows results for the S+S reaction at 200 GeV/nucleon bombard- 
ing energy together with the experimental data. Column 1 displays multiplic- 
ities observed in experiment NA35 |]l6l, fl^. Column 2 contains the hadron 
multiphcities predicted by the ALCOR model. Adjusting the three above pa- 
rameters of ALCOR one obtains fairly good agreement with the experimental 
data for particles ir' , , K~ , K^,Y^ and y". (Here we introduced the no- 
tation F° = A° + S° and F° = a" + if for the measured neutral strange 
and anti-strange baryons .) We took into account that the number of partici- 
pant nucleons was measured to be N^^^^-^ = 51 |TB[. We used the parameters 
^tot,pair = 158.1 with gs = 0.255 for the strangeness production, and gs = 0.04 
for the baryon formation . For comparison we display some results of the HI- 



JING, RQMD and QGSM models in Columns 3, 4 and 5 taken from refs. |T4 
[0 (early non-collective string model versions). Table 2 shows results for the 
strange baryon and anti-baryon ratios. Column 1 displays the experimental 
data [0 and Column 2 contains the ALCOR predictions. The agreement 
seems to be good. 



S+S 


WA94 


ALCOR 


yo/yo 


0.23 ±0.01 


0.23 




0.55 ±0.07 


0.50 




0.09 ±0.01 


0.11 




0.21 ±0.02 


0.24 



Table 2: Strange baryon and anti-baryon ratios measured by WA94 Collab- 



oration 1 19] and obtained from ALCOR for S + S collision at 200 GeV/nucl 



bombarding energy. 

From this experience we conjecture that the above rehadronization process 
is locally quick: the quarks and antiquarks (including the fragmented gluons) 
appear in hadrons according to simple kinematic rules and production branch- 
ing ratios. This finding, however, does not settle the question whether quark 
matter with any degree of collectivity has been formed or other hadronization 
processes including strings or color ropes are the source of this amount 
of quarks. Local thermal distribution of the different quark flavors and the 
presence of a longitudinal flow cannot, on the other hand, be excluded on the 
basis of this experimental data. 
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To obtain results for other reactions, we assume a scaling for the produced 
quark-antiquark pairs as 

Nttani^) = KUr. (v^) (17) 

\ ^partic / 

where NpJ:^^^ are the number of participant nucleons in the A+A collision. 

The scaling exponent a may have the value a = 1, or, for more collective 
production processes one expects a = 4/3. Furthermore we shall use the values 
for qb and Qs obtained above. If one of these parameters has to be changed in 
order to approve the correspondence of the calculated particle numbers with 
the experimental ones, than this change must have a physical interpretation. 
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0.013 ±0.004 


0.019 



Table 3: Hadron multiplicities for p + n and p + p collision observed ex- 
perimentally [^ and obtained in the ALCOR models at 200 GeV/nucleon 
bombarding energy. The numbers in square brackets are the average of the 
corresponding p + p and n + n results. The numbers in normal brackets are 
assumed values with 50% uncertainty [EO. 
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With these assumption we calculated the particle production also for nuc- 
leon-nucleon collisions namely for proton-proton and proton-neutron collisions 
at 200 GeV/nucleon energy. For this case, in order to get the number of 
strange hadrons correctly, we had to decrease the strangeness production fac- 
tor. Table 3. shows the ALCOR results for p + n and p + p collisions to- 
gether with the experimental values from Ref. In both cases we used 
the same parametrization, namely N^J^^^^ = 2, N^l^^^^^ = 4.5 with a decreased 
strangeness formation factor Qs = 0.16 and an unchanged baryon formation 
factor gB = 0.04. The scahng exponent is a = 1.1. (In the two cases only 
the flavour component of the incoming participant nucleons were different, i.e. 
p + p and p + n, respectively.) One can observe, that the calculated meson 
and non-strange baryon numbers are not far from the measured ones. To 
achive this we had to decrease the strangeness production factor, gg, from its 
value in S* -|- S* reaction. Thus we may conclude, that in the S + S reaction 
the strangeness production is enhanced with respect to the nucleon- 
nucleon reaction by a factor of 1.6. Furthermore we may observe from 
Table 3., that the calculated y° = S° -|- A'' and = + A*^ numbers deviate 
essentially from the measured ones. Namely ALCOR overpredicts the produc- 
tion of these particles. This discrepancy can be interpreted by assuming that 
the hadronization mechanism in the S + S reaction is different from 
that for the nucIeon-|-nucIeon reaction and this can be seen undoubt- 
edly on strange baryon production. In fact, the ALCOR model assumes, that 
there may be quark exchanges between neighbouring excited objects (strongly 
packed strings, etc.) allowing larger role to the combinatoric possibilities. 
Thus one may describe this formation of a semi-deconfined state 
of the matter. However, in nucleon-nucleon collision this effect is missing. 

Using the parametrization of S-l-S collision and the above scaling law, we 
may make predictions for the Pb-|-Pb collision at 160 GeV/nucleon energy. In 
this case the number of participant nucleons is A^^^f* = 390 ± 10 obtained 
from Monte-Carlo simulations [^. We will use the mean value. During the 
extrapolation of the total number of produced quark-antiquark pairs from the 
S-l-S collision one needs to consider an energy rescaling. We will assume a 
logaritmic one: 

tot,pair V 



This rescaling yields a ~ 4% correction and one obtains N^Jl^^^^ = 1164 in the 
linear case (a = 1) and N^^^^-^ = 2294 in the more collective one (a = 4/3). 
Table 4 displays the predictions of ALCOR for both scaling cases and the last 
column contains some results from HIJING taken from ref. ill . 
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Table 4: Hadron multiplicities predicted by the ALCOR model assuming 
different scaling laws discussed in the text. For comparison, predictions of the 
HIJING model taken from ref. 



141 are also selected. 
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Conclusion 



We have applied the hadronization model, ALCOR, for the analysis of the 
experimental data of p + p, p + n and S + S collisions and for the prediction 
of hadron multiplicities in the Ph + Ph collision at CERN SPS energy. 

This model uses dynamical coalescence probabilities obtained from hadro- 
production reaction rate factors obtained on a one dimensional scaling flow 
background. The model is sensitive to only three free parameters: the total 
number of quark-antiquark pairs produced, an in-mcdium baryon formation 
factor and the strangeness production factor (7^, fitted to the experimental 
results of the reaction S+S. (The efi:ect of changes in common factors are 
removed by the normalization constants.) With no further adjustment we have 
been able to describe several observed hadron multiplicities in agreement with 
the NA35 and WA94 experiments. We could clearly identify the strangeness 
enhancement m. S -\- S collision with respect to the nucleon-nucleon reaction. 

Furthermore, for p + p and p + n collisions the ALCOR model could not 
reproduce the measured strange baryon and antibaryon numbers on the ba- 
sis oi S -\- S coUision. We interpret this discrepancy as a signature showing 
that the hadronization mechanism in the S + S reaction differs from that for 
the nucleon+nucleon reaction. Wc infer from this comparison that in S + S 
collision a sort of semi-dcconfincd state of the matter was formed which can 
be characterized by enhanced quark exchanges between neighbouring excited 
objects. This is a basic assumption in the ALCOR model. 

We made detailed predictions for hadron multiplicities in a Pb+Pb reaction 
assuming two different scaling exponent. There is an agreement between the 
HIJING and the ALCOR model (in hnear case, a = 1) for prediction of total 
negative particle, pion and kaon multiplicities. However, the strange baryon 
and anti-baryon yields differ from the HlJlNG result. Now we look forward to 
compare these predictions with the results of the NA49 CERN SPS experiment 
to determine the details of hadronization processes. Possible agreement would 
support our statement about the formation of semi-deconfined state of the 
matter in high energy heavy ion collisions. 
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